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"Naked" Cronin effect in A+A collisions from SPS to RHIC 
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Abstract Baseline computations of the Cronin effect in nuclear collisions at energies spanning the SPS and 
the RHIC accelerators are performed in the Glauber-Eikonal model, which ascribes the effect to initial- 
state incoherent multiple parton scatterings. The model accounts very well for the mid-rapidity Cronin 
effect in hadron-nucleus collisions in the y/s — 27 — 200 GeV center of mass energy range, and will be 
extended to nucleus-nucleus collisions. The computations are performed under the assumption that the 
partons do not interact with the medium produced in the collision. Therefore, medium effects such as 
energy loss in a Quark-Gluon Plasma may be detected and measured as deviations from the presented 
baseline computation of the "naked" Cronin effect. 

PACS. 12.38.Mh - 24.85. +p - 25.75.-q 



1 Introduction 

The Cronin effect £Q is the deformation of hadron trans- 
verse momentum spectrum in hadron-nucleus (p+A) and 
nucleus-nucleus (A+A) collisions relative to linear extrap- 
olation from proton-proton (p+p) reactions. The Cronin 
effect may be quantified by taking the hadron transverse 
momentum (pr) spectrum in p+A collision in a given cen- 
trality class (cel.), normalizing it to binary scaled p+p 
collisions by the inverse thickness function Ta, and finally 
dividing it by the p+p spectrum: 



RhA = 



1 



Ta{c.cI.) dp\,dy 



dN pA-,hX 

(cel.) 



da pp - 
dp\dy 



(1) 



For A+A collisions, an analogous definition of Raa is 
used. To cancel systematic errors as much as possible, it 
is also customary to take the ratio of a given centrality 
class to the most peripheral one. 

In p+A collisions, the ratio RhA is smaller than 1 at 
low momentum pt^2 GeV (the spectrum is suppressed). 
At intermediate pt^2-6 GeV is larger than one (the spec- 
trum is enhanced) and finally tends to 1 at higher pr's 
(the spectrum follows a p+p scaling). The described bell 
shape of the "Cronin ratio" RhA has been observed on a 
wide energy range s/s sw 20 — 200 and for different hadron 
species PE]. 

In A+A at low energy ^ = 17 A GeV at SPS 0, 
one observes a similar behaviour of the Cronin ratio as in 
p+A collisions. However, at RHIC energy y/s = 63 — 200 
GeV, the above described bell shape appears on top of a 
largely suppressed spectrum 0]. This suppression is often 
interpreted as due to energy loss of a parton traveling 
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a medium of density up to 15 times the normal nuclear 
density j^j, largely sufficient for the creation of a Quark- 
Gluon Plasma (QGP). 

An intriguing fact in both the p+A and A+A cases 
is the much larger magnitude of the effect on baryons 
compared to mesons, with the difference between the two 
seeming to decrease with energy P]E]- This aspect of the 
famed "baryon/meson" anomaly in high energy nuclear 
collision has not yet found an accepted explanation in the- 
ory, and will later come back into our discussion. 

In this note, I discuss the Glauber-Eikonal (GE) model 
computation [2] of the Cronin effect in A+A collisions. 
The model describes multiple parton scatterings in pQCD 
and extrapolates the physics of proton-proton (p+p) col- 
lisions to p+A and A+A collisions without ad hoc param- 
eters for the nuclear case. This allows a baseline compu- 
tation of the Cronin effect in A+A without the further 
complications of medium induced modifications of the pr 
spectrum. This is what I call the bare, or "naked" , Cronin 
effect. Comparison with experimental data allows to mea- 
sure the medium induced effects, and in particular to de- 
tect when the medium begins to suppress hadron spectra 
as a function of the energy and centrality of the collisions. 



2 The Glauber-Eikonal model 

The Glauber-Eikonal (GE) approach [Jj to the Cronin ef- 
fect treats multiple 2— >2 partonic collisions in collinearly 
factorized pQCD. The low-px spectra in nuclear collisions 
are suppressed by unitarity. At moderate pr, the accumu- 
lation of transverse momentum leads to an enhancement 
of transverse spectra. At high pt the binary scaled p+p 
spectrum is recovered. 

The cross-section for the production of a hadron h with 
transverse momentum pt and rapidity y in p+A collisions 



2 



Alberto Accardi: "Naked" Cronin effect in A+A collisions from SPS to RHIC 



at fixed impact parameter b is written as 



da 



pA-^hX 



cPpxdycPb 



= (xfi/p) 



da 



iA 



T A (b)(xf j/A )(g> 



d 2 pTdyid?b 
da jp 



D 



(2) 



d 2 p T dy 3 



D 



where the crossed-circle symbols denote appropriate inte- 
grations and summations over parton flavours i and j, see 
Ref. [3| for details.. The first term in Eq. (j2J accounts for 
multiple semihard scatterings of the parton i on the target 
nucleus; in the second term the nucleus partons are as- 
sumed to undergo a single scattering on the proton. T A {b) 
is the target nucleus thickness function, and is the 

fragmentation function of the i parton. The average par- 
ton flux from the proton, (xfi/ p ), and the parton-nucleon 
cross-section, da lN , are defined as 



da lN da i] 

da 1 



(3) 



(xfi/p) = K ii/i/p ® —z (po) ® Xjf j/N ((k^)) 



( da 



\d 2 p T dyi 



(4) 



where only the sum over j is understood, t is the Mandel- 
stam variable and da/dt are leading order parton-parton 
cross-sections in collinearly factorized pQCD. K is a con- 
stant factor which takes into account next-to-leading order 
corrections. To regularize the IR divergences of the single- 
scattering pQCD parton-nucleon cross-sections, a small 
mass regulator pq is introduced in the propagators, and 
Q = y/Pr + Po/2 is the scale of the hard process. Finally, 
a small intrinsic transverse momentum (fef,) = 0.52 GeV 2 
is introduced to better describe the hadron spectra in p+p 
collisions at intermediate pr ~ 2 — 5 GeV. Since all pa- 
rameters have been fixed in p+p collisions, we are able to 
compute the spectra in p+A and A+A collision with no 
extra freedom. 

The free parameters po and K in Eqs. are fitted 

to, viz., low- and high-pT hadron production data in p+p 
collisions at the energy and rapidity of interest. For this 
study, I took from Ref. the values of po extracted at 
y/s = 27 A and 200 GeV, viz., p Q = 0.7 and 1.0 GeV, and 
assumed a logarithmic dependence on the energy: 



po = 0.151 + 0.100 log f— ^— ) ±10% GeV. 

V 1 GeV / 



(5) 



As regards the K factor, it turns out to be K sa 1 quite 
insensitive to the energy of the collision in the 27-200 GeV 
energy range, and I will use this value throughout this 
note. Clearly, one would need more systematic studies of 
spectra in p+p collisions, especially when extrapolating 
the parameters at the SPS energy of 17 GeV. These will 
be presented elsewhere. 

Nuclear effects are included in da lA , the average trans- 
verse momentum distribution of a parton who suffered at 





p. (GeWc) 



Figure 1. Cronin effect on neutral pion production as a func- 
tion of centrality in d+Au collisions at RHIC. Curved lines 
and band are the GE model computation with its theoreti- 
cal uncertainty. The horizontal gray band is the experimental 
uncertainty in the normalization of the ratio. Figure adapted 
from Ref. 0. 



least one semihard scattering on the target nucleus A: 



da 



iA 



d 2 pTdyd 2 b 



jrl. j d 2 bd 2 k 1 ---d 2 k n s( ki- PT ) 



71=1 

da m 
d 2 ki 



x e 



T A (b) x ■ • • x 



r (po)T A (b) 



da 



iN 



d 2 k, 



-T A {b) 



(G) 



This equation sums all processes with n multiple 2— >2 par- 
ton scatterings. The exponential factor in Eq. © repre- 
sents the probability that the parton suffered no semihard 
scatterings after the n-th one, and explicitly unitarizes the 
cross-section at the nuclear level. 

The above presented GE model describes quite well the 
energy and centrality dependence of the Cronin effect on 
pion production in p+A collisions at Fermilab y/s — 27 A 
GeV, and in d+Au collisions at RHIC ^fs = 200 GeV, see 
Ref. 7 and Fig.^ As all other models based on multiple 
initial-state parton scatterings, the model fails in describ- 
ing the effect on baryon production. The main reason is 
that baryon production is not understood in this frame- 
work even at the p+p level, which is the starting point 
of the GE approach, see Eqs. J3J and J3J. An alterna- 
tive explanation of the baryon/meson anomaly based on 
final-state parton recombination |10| is more successful in 
this respect than the independent fragmentation approach 
used in the GE model. 



3 The "naked" Cronin effect in A+A 
collisions 

The natural generalization of Eq. J2J to the case of nucleus- 
nucleus collisions is to allow the partons from the target 
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Figure 2. Naked Cronin effect on ir° production in Au+Au 
collisions at ^fs — 62 GeV compared to preliminary PHENIX 
experimental data from 0. 



nucleus to rescatter on the projectile nucleus, as well: 



da 



AB^hX 



= d 



cPpxdycPb 

x S^T A (b-s){xf i/A ) 

+ T B (s)(xfj/ B ) 



(7) 



da 



iB 



d 2 pTdyid 2 s 
da lA 



d 2 p T dy l d 2 (b - s) 



D 



Vh=-y . 



As no medium effects arc included in this formula, what 
we are computing is the "naked" Cronin effect, i.e., the 
effect stripped of medium-induced modifications such as, 
e.g., jet quenching. Hence, the resulting Cronin ratio may 
be used as a baseline to detect and measure jet quenching. 

Consider, e.g., the preliminary PHENIX data on ir° 
production in Au+Au collisions at y/s = 63 GeV in Fig.0 
The characteristic bell shape of the Cronin effect at mod- 
erate pt is clearly suppressed. However, the magnitude 
of the suppression is not self-evident from experimental 
data alone: a computation is needed of the unsuppressed 
spectrum. The comparison of experimental data and the 
GE model computation shows a quite large jet quenching. 
Going from central to peripheral collisions, the suppressed 
experimental Cronin peak and the unsuppressed theoreti- 
cal curve approach each other. As we can expect this trend 
to continue in more and more peripheral events, a natu- 
ral question arises: at what centrality does the quenching 
starts? This is clearly a very important question, whose 
answer will help to pinpoint the conditions for QGP for- 
mation. 

As PHENIX data on ir° in more peripheral bins are not 
yet available, we may turn to PHOBOS data on charged 
hadron (/i ± ) production. Here we come across the bar- 
yon/meson anomaly previously discussed: the GE model 
underestimates the Cronin effect on protons and anti- 
protons, which give a sizable contribution to the h^ 1 sam- 
ple. To attempt nonetheless an answer to the question of 
when does the medium start to quench hadrons, we may 
compare data to the computation for 7r°'s, see Fig- EJ 
The required modification of the shape and magnitude 
of theoretical curves can be estimated by comparing the 
PHOBOS data points with PHENIX tt° data point in 



similar centrality bins. With the above caveats in mind, we 
can estimate that hadron quenching starts at a centrality 
of around 40%. 

I would like to stress the importance of the availability 
of a theoretical computation of the naked Cronin effect on 
the example of the 4 most peripheral bins in Fig. (up to 
15-25% centrality). The data points do not show a medium 
effect by themselves: they arise above 1 at pt ~ 1.5 — 2 
GeV, display a nice Cronin peak, and seem to converge 
to 1 again at higher px inside the systematic errors. This 
behaviour would point to an unsuppressed Cronin effect, 
but the comparison to the naked Cronin curve (modulo 
the baryon/meson anomaly) shows a sizable medium ef- 
fect starting at around 40%, which otherwise might have 
passed unnoticed. 

Finally, we can address the WA98 data on ir° produc- 
tion in Au+Au collision at y/s = 17 GeV [3J. Here we 
have two complications. On the experimental side, there 
is a large uncertainty on the p+p spectrum to be used as 
a normalization in Eq. (TJJ) , due to different extrapolations 
of higher energy data, see Ref. for more details and 
references. On the theory side, we have an uncertainty on 
the poi^fs) and K(tJs) parameters due to the extrapola- 
tion at this low energy of values extracted [7] from data at 
\fs = 27 and 200 GeV. Preliminary computations show a 
moderate quenching in central collisions. 



4 Perspectives 

A cross check of WA98 results may be obtained using the 
minimum bias data sample of NA60 ^2 for In+In colli- 
sions at y/s = 17 GeV. The dimuon trigger data, which 
has a bias toward central collisions, may also be used. 

NA60 can also perform a nice experimental study of 
the system size dependence of the Cronin effect on ft+ 
production in p+A collisions at low y/s = 17 — 27 GeV. 
Indeed, the experimental setup has 7 nuclei (ranging from 
beryllium to uranium) mounted in the target region and 
sharing the same beam ^21 • This allows to measure the ra- 
tio R a of differential cross-sections on two different nuclei 
because the beam luminosity cancels in the ratio 



r: 



pA/pB 



1 da pA ~ 
A dp\dy 



■ i,x 



1 da pB ^ hx 



B dpipdy 



(8) 



where B is the lightest available nucleus, in this case beryl- 
lium. The A dependence of the ratio is analogous to the 
centrality dependence, but eliminates the large experimen- 
tal uncertainties due to the determination of the centrality 
and to the normalization of the Cronin ratio. Moreover, 
without need of centrality cuts the statistics may be suffi- 
cient to probe the high-pT region and test the multi parton 
scattering mechanism, which predicts a divergence of the 
Cronin ratio at pr — 4 GeV and 8 GeV in collision of, viz., 
y/s = 17 GeV and 27 GeV. Theoretical computations will 
need to take into account the bias toward central collisions 
induced by the dimuon trigger, which is the only trigger 
used in p+A collisions. 
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The HERMES data on nuclear DIS can shed light on 
the baryon/meson anomaly and the space-time evolution 
of the hadronization process in a cleaner and more control- 
lable setting than nuclear collisions. I would like to invite 
heavy ion theorists and experimentalists to carefully read 
HERMES papers: there is a lot to be learned which is di- 
rectly relevant to RHIC physics. 



The baryon/meson anomaly is a pervasive phenomenon 
in QCD: it is not limited to hadronic and nuclear colli- 
sions, but present also in the HERMES measurements of 
hadron production in deep inelastic scattering on nuclei 
(nDIS) ^21 • I n this case, the ratio of hadron multiplicities 
on a nuclear target and on deuterium is measured as a 
function of z — Eh/v, where Eh is the hadron momen- 
tum and v is the virtual photon energy. Experimentally, 
meson quenching is observed on the whole z — 0.2 — 1 
range, while protons are enhanced at z < 0.4 and sup- 
pressed above; both the quenching and enhancement in- 
crease with A. At HERMES the typical hadron momen- 
tum Eh « 2—12 GeV is comparable to midrapidity hadron 
momenta measured at RHIC: there is a lot to be learned 
from nDIS about hadronization, both experimentally and 
theoretically, which can be applied to RHIC physics! 

5 Conclusions 

The Glauber-Eikonal model is successful in explaining the 
Cronin effect on midrapidity meson production in p+A 
collisions on a wide range of energy, ^/s =27-200 GeV. 
When extended to A+A collisions, it gives a reliable base- 
line computation of what the Cronin effect should look 
like in absence of the medium-induced hadron quenching. 

By comparing the GE computation with experimen- 
tal data it is possible to study the onset and magnitude 
of hadron quenching. As an example, I studied recent 
PHENIX and PHOBOS data in Au+Au collisions at y/s = 
63 GeV, and found that the produced medium begins to 
quench hadron spectra in collisions of around 40% cen- 
trality, with a strength quickly increasing with centrality. 
Preliminary analysis of WA98 data at y/s = 17 GeV sug- 
gests a moderate hadron quenching in central collisions. A 
systematic analysis of all the available experimental data 
is under preparation. 

More can be learned about the Cronin effect in low- 
energy p+A collisions from the unique multiple target 
setup of the NA60 experiment by studying, as proposed in 
this note, the A-dependence of the ratio (jHJ of differential 
cross-sections on a target nucleus A and on beryllium, the 
lightest available target. 
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